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We demonstrate that a quantum spin Hall current, spontaneously generated at the edge of a two- 
dimensional topological insulator, acts as a source of spin-pumping for a magnetic impurity with 
uniaxial anisotropy. One can then manipulate the impurity spin direction by means of an electrical 
current without using either magnetic electrodes or an external magnetic field. Furthermore we 
show that the unique properties of the quantum spin Hall topological state have profound effects on 
the inelastic electron tunneling spectrum of the impurity. For low current densities inelastic spin-flip 
events do not contribute to the conductance. As a consequence the conductance steps, normally 
appearing at voltages corresponding to the spin excitations, are completely suppressed. In contrast 
an intense current leads to spin pumping and generates a transverse component of the impurity 
spin. This breaks the topological phase yielding to the conductance steps. 



The Quantum Spin Hall (QSH) state is a new topolog- 
ical phase of matter, which has attracted growing atten- 
tion over the past few years. It is time-reversal invariant 
and characterized by both a bulk gap and gapless helical 
edge states with opposite spins counter-propagating at a 
given edge. In a pathfinding work by Kane and Mele p] 
it was shown that at low energy a QSH effect can be gen- 
erated in a sheet of graphene subject to spin-orbit inter- 
action. A topological invariant, Z2, distinguishes such a 
QSH insulator from an ordinary one [2]. Unfortunately, 
spin-orbit interaction in graphene is too small to prac- 
tically realize the QSH state at a realistic temperature, 
however a few other options are available. Bernevig et al. 
predicted that HgTe/CdTe quantum wells might exhibit 
this novel phase [3], a prediction soon confirmed by ex- 
periments [4, 5j . Three-dimensional analogs of QSH state 
have also been found and are generically termed topolog- 
ical insulators (TIs). More recently, evidence for helical 
edge modes in InAs/GaSb quantum wells has also been 
found [5J. Interestingly, both Silicene and two dimen- 
sional Ge have been predicted to exhibit QSH state at 
experimentally accessible temperatures [TJ [8] , and there 
are theoretical proposals to realize the same in graphene 
by using non-magnetic adatoms [5]. These last two op- 
tions would yield a material with a low-energy electronic 
structure identical to that proposed by Kane and Mele. 

Given the peculiar spin structure of the QSH state it 
becomes natural to ask ourselves whether this can be 
used to manipulate magnetic objects [ID] . In particular 
the question we address in this letter is whether spin- 
pumping at the single spin level can be achieved without 
using spin-polarized electrodes or an external magnetic 
field. In a nutshell we wish to propose an analog to 
the numerous recent investigations concerning spin-flip 
inelastic electron tunneling spectroscopy (SF-IETS) for 
magnetic adatoms on insulating surfaces [11| . either in 
equilibrium or in spin-pumping conditions |12j . 

Here we demonstrate that a magnetic impurity de- 
posited at the edge of a Z2 TI and presenting a uniaxial 
magnetic anisotropy, which makes it non-Kondo-active 



|13) . can be manipulated by the QSH current. Further- 
more we show that the topological nature of the QSH 
state has profound consequences on the SF-IETS con- 
ductance spectrum. At low current intensity there is a 
complete suppression of the conductance steps appearing 
at the critical biases characteristic of the activation of an 
inelastic spin-scattering channel In contrast, for cur- 
rents large enough to produce spin-pumping the spin of 
the magnetic impurity is driven away from the anisotropy 
axis. This breaks the topological protection of the helical 
edge states and the conductance steps reappear. Intrigu- 
ingly, in this situation there is a strong dependence of the 
SF-IETS conductance spectrum on the bias polarity. Our 
calculations are conducted by using the non-equilibrium 
Green's function method for transport combined with a 
perturbative approach to spin-scattering from magnetic 
impurities [MHTB] , 

The device that we have in mind is schematically pre- 
sented in Fig. [T] and consists of two semi-infinite cur- 
rent/voltage electrodes sandwiching a Z2 TI ribbon in 
which a magnetic impurity is absorbed at one of the two 
edges. Our working hypothesis is that one can construct 
such a device with either a strong or a weak electronic 
coupling between the electrodes and the ribbon, i.e. that 
the interface resistance can be engineered. The entire 
system is described at the tight-binding level and for the 
electrodes we use a simple square lattice with hopping 
parameter, t v , (77 = L, R). 

The TI ribbon has a honeycomb lattice structure with 
zig-zag edge geometry and it is described by the 2D Kane- 
Mele (KM) Hamiltonian, -ffxivb which reads 

HkM = e ^ £,ic\ a C ia + t\ ^ 4a C ja+ (1) 
ia {ij) , a 

+ it 2 ^ Vij^WUpCja ■ 
{(ij)), a/3 

The first term describes a staggered sublattice potential 
with on-site energy £o and & being £j = +1 for the A 
sublattice and — —1 for the B one. Here c\ a (c ia ) 
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FIG. 1: Schematic representation of the device considered in this 
work (a), comprising a TI with honeycomb lattice structure and 
a magnetic impurity adsorbed at one of the two edges. In the 
transport calculations the electronic structures of the electrodes is 
described by the self-energies, £„ (n = L, R). In panel (b) the car- 
toons describe a spin-flip scattering event. A right-going electron 
with up spin direction (left panel) is inelastically back-scattered 
by the magnetic impurity. In the process both the electron and 
the impurity spins are reversed (right panel). Note that, given the 
topological nature of the ribbon, spin-flip forbids electron trans- 
mission as the edge presenting a right-going spin-up state does not 
possess a right-going spin-down one. 



creates (annihilates) an electron at site i with spin a. 
The second term is the nearest neighbour hopping with 
strength t\ {t\ sets the energy scale of the problem). Fi- 
nally the third term, which drives the topological phase, 
is a second nearest neighbour hopping with strength £2 
(i = v— 1) . This describes the coupling of the electrons 
orbital motion to their spins via the z-component of the 
Pauli matrices (<r z ). The parameter Vij is +1 for counter- 
clockwise hopping and —1 for clockwise. 

If we now attach an impurity at site /, the total elec- 
tronic Hamiltonian will become 



H 



KM 



E 



c \ a c ia i (2) 



where in addition to Hkm one has the on-site potential of 
the impurity, £/, and the hopping between the impurity 
site / and its neighbor i on the honeycomb lattice (with 
strength tj). Finally there are two other terms related to 
the magnetic impurity spin, S 

H sp — DSl ; H e l- a p = J s d c I a [ tT }al3 c iP ' ^ • (3) 

a/3 

The first, H sp , describes the uniaxial anisotropy (along 
z) with D being the zero-field splitting parameter. The 
second, H c i_ sp , couples the current-carrying electrons to 
the local impurity spin with interaction strength J sc j (cr 
is a vector of Pauli matrices). This is usually known as 
the s-d model for magnetism |17j . 



Electron transport is investigated within the 
Landauer-Biittiker approach |18j implemented by the 
non-equilibrium Green's function (NEGF) method |19j . 
The central quantity to evaluate is the retarded 
electronic Green's function for the scattering re- 
gion (the TI ribbon) in presence of the electrodes, 
G r = [{E + iO+)I - H cl - E L - Sr.]- 1 , where 
(r) = L, R) are the electrodes self-energy, which can be 
computed with standard techniques [30]- These depend 
of the hopping parameter between the ribbon and the 
electrodes, t v R,, whose magnitude sets the intensity of 
the current. 

When the conducting electrons couple to the impurity 
spin (J sc j ^ 0), the problem becomes intrinsically many- 
body in nature. This is made treatable by constructing 
a perturbation theory in the electron-spin Hamiltonian, 
which allows us to incorporate the effect of the electron- 
spin interaction in an additional self-energy, £i n t. In 
this work we truncate the perturbation expansion to the 
second order [TU [TB] in both the electron and the im- 
purity spin propagator. The latter contains informa- 
tion about the state of the magnetic impurity spins, 
through the population, P n , of the eigenvectors of the 
spin-Hamiltonian, H sp . In particular it is possible to 
show that the P ra 's satisfy a master-equation of the form 



E 



[P n (l-Ih)W ln -I\0.-P n )W n t]+(I%-P n )/P, 



dP n 

dt 

(4) 

where the detailed expression for the transition rates, 
Wi n7 can be found in Ref. [IB] , and j3 = k^T with 
being the Boltzmann constant and T the temperature. 
In Eq. [i] the populations P® are those of the ground 
state. With this at hand we can compute the current, 
I, and hence by numerical derivative the conductance, 
G = dl/dV. 

We start our discussion by first looking at the trans- 
port properties of the ribbon in absence of the magnetic 
impurity. The relevant quantity here is the spin-resolved 
total transmission coefficient along a particular edge [21] , 
which is given by 



^(e F )=Tr„ x [r L G"T R G a ]^ 



(5) 



where a is the spin index (a =1~>4-)) s labels the edges 
(s = top, bottom) and G a is the advanced Green's func- 
tion. The trace is over the number of atoms, n x , along the 
given edge and the transmission coefficient is evaluated at 
the Fermi energy, ep. As a matter of notation a {n x ,n y ) 
ribbon contains n x atoms in the direction of transport 
and n y along the transverse one. When the Fermi level 
is fixed at the half-filling point the ribbon is insulating in 
the bulk, but presents edge topological protected states 
(here = t v = t\ = 1, and t% = ii/3, which ensures 
that the KM Hamiltonian describes a QSH state). In 
this situation we find for a (11, 6) ribbon, T*° p = 0.9, 
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T™ v = 0.1, T^° ttoln = 0.1 and T°° ttom = 0.9. Such val- 
ues demonstrate that the current along the QSH edges 
is spin-polarized, although not completely because of the 
finite size of the ribbon. Calculations for a (7, 4) ribbon 
give us T*° p = 0.85, T*J P = 0.15, T t b ° ttom = 0.15 and 

T bottom = o g5 

We now switch on the magnetic interaction between 
a S = 1 local spin and a (11, 6) ribbon. In general we 
place the impurity at the centre of the edge and choose a 
coupling parameter, tj, and an onsite energy, si, to en- 
sure that the density of states localized at the impurity 
site, pi(E), is approximately constant for energies, E, 
around the Fermi level (this ensures the convergence of 
the perturbation scheme [2]). The exchange coupling, 
J s d, is chosen so that the perturbation parameter, pj J s( j, 
is approximately 0.1. These conditions are satisfied for: 
cj = J s d = h/2 and tj — ii/4. The spin degeneracy 
is lifted by an axial anisotropy D = — 10~ 3 t\, which 
corresponds to D = —2.0 meV, if the nearest neighbour- 
hopping in the ribbon is fixed at a reasonable value of 
<i = 2 eV (fceT = 0.05). The uniaxial anisotropy gives 
us a degenerate ground state with the two spin states 
| — 1) and | + 1) separated from the first excited state 
|0) by \D\. As a result we do not expect a Kondo-like 
behaviour since no allowed transition between the de- 
generate ground state may occur. The second order per- 
turbation expansion is then well justified. The values 
trj = 4i, ; R = t\ ensure that the spin system remains in 
equilibrium, i.e. in its ground state, throughout the spin- 
scattering process. 

Figure [2] shows the calculated conductance spectra, 
G(V), normalised to the V = conductance, G e i, for 
three values of the parameter governing the QSH state, 
For t2 — there are no topologically protected edges 
and we observe the characteristic inelastic conductance 
step at a voltage V = D/e, when the transition from the 
ground state to |0) becomes possible (e is the electron 
charge). However, as is increased and we enter into the 
topological phase, we reveal a suppression of the inelastic 
contribution to the conductance, with an almost full sup- 
pression at the maximum value of t<x = ti/3. The cartoon 
in Fig. |T|(b) helps to understand the mechanism for such 
a suppression. At a positive bias, the right-going current 
is up spin-polarized. This means that the | — 1) — > |0) 
transition scatters out spin-down electrons. These can- 
not propagate towards the right electrode since there is 
no right-moving spin-down state in the upper edge and, 
as a consequence, they are completely reflected. Hence, 
as spin-flip events can only lead to backscattered elec- 
trons, the inelastic channel does not contribute to the 
current. Note that the residual conductance increase in 
Fig. [2] for £2 = ti/3 is simply due to the finite size of the 
ribbon. 

We now investigate the possibility of manipulating the 
impurity spin direction. This is achieved by increasing 
the overall conductance, i.e. by increasing the average 
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FIG. 2: (Color online) SP-IETS conductance spectrum for a TI 
(11, 6) ribbon with a S = 1 magnetic impurity attached at the up- 
per edge. Note that the conductance step at the voltage character- 
istic of the inelastic excitation gets suppressed as the t2 parameter 
is increased, i.e. as the ribbon is brought well inside the topological 
region of the phase diagram. 

current density. When one works with an STM setup 
bringing the tip closer to the impurity [TJ] does the job, 
while here our control parameter is the electronic cou- 
pling between the leads and the ribbon, t V R. As such 
all the calculations that follow have been performed with 
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FIG. 3: (Color online) (a) Non-equilibrium population as a func- 
tion of bias of the 5=1 impurity spin states for a (7, 4) (broken 
lines) and a (11, 6) ribbon (solid lines). In panel (b) we show the 
average magnetization of the impurity for the same ribbons. 

The calculated populations of the various spin states 
are plotted as a function of bias in Fig. (3^a) for both 
a (11, 6) and a (7, 4) ribbon. AS=1 spin in a uni- 
axial anisotropy field and in thermal equilibrium with 
an electron bath presents an equal probability to oc- 
cupy the I + 1) and the | — 1) states, i.e. for V = 
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FIG. 4: SP-IETS conductance spectrum for a TI (11, 6) ribbon 
with a S = 1 magnetic impurity attached at the upper edge. In 
this case the current is intense and drives the impurity spin away 
from the uniaxial anisotropy axis (see Fig.[3]|. Notably now there is 
a step in the differential conductance at the voltage corresponding 
to the inelastic transition | ± 1} — > |0), The magnitude and sign of 
such step depends on the bias polarity. In the inset the inelastic 
contribution to the conductance. 



one has P + i = P_i = 1/2. As soon as the bias is in- 
creased at and above \D\/e, excitations to the |0) state 
become possible due to spin-flip back-scattering. In this 
case however the current is intense, so that in between 
two scattering events the impurity spin does not have the 
time to relax back to the degenerate ground state. This 
means that now a spin-up electron (the majority specie 
in the upper edge right-going channel) can also induce 
the transition |0) — > \ + 1). The consequence is that the 
electronic current flowing at the upper edge, in virtue of 
its spin polarization and its intensity, produces a net flow 
of population between the two degenerate ground state, 
i.e. for V > +\D\/e one has P+i > P_i. In other words 
the impurity spin is driven by the current away from its 
uniaxial anisotropy axis. This can be fully appreciated 
by looking at Fig[3]jb), where we show the average mag- 
netization (S z ) = J2 m P m S* n as a function of bias. Such 
spin-pumping is essentially identical to what happens for 
spin-polarized tips |22| except that now one does not need 
either a magnetic electrode or an external magnetic field. 
Note that at a negative bias the effect is reversed, i.e. for 
V < —\D\/e one has P_i > P+i, and that placing the 
impurity on the lower edge is equivalent to reversing the 
bias polarity. 

The effects of the spin-pumping on the shape of the 
conductance spectrum are finally presented in Fig. [4] 
This time the G(V) trace presents a step at the voltage 
corresponding to the | ± 1) — > |0) transition, i.e. the elec- 
tron transport becomes sensitive to spin-flipping events. 
Such an appearance of the conductance step signals the 
suppression of the topological helical states induced by 
the transverse magnetization of the spin impurity |23j . 



Intriguingly, the magnitude and sign of the conductance 
step depends on the bias polarity. In particular we note 
that there is an inelastic contribution, which is symmetric 
with respect to the sign of V, and an elastic one which 
is anti-symmetric, i.e. the elastic current increases for 
V > \D\/e and decreases for V < — \D\/e. Placing the 
impurity on the opposite edge yields a mirror symmetric 
conductance trace (i.e. V — > —V). 

In conclusion we have demonstrated that a QSH cur- 
rent flowing at the edge of a TI can be used to ma- 
nipulate the spin of a magnetic impurity. This does not 
require either an external magnetic field or magnetic elec- 
trodes, i.e. it allows one to implement spintronics with- 
out magnetism. Importantly the fingerprint of the ma- 
nipulation can be found in the conductance profile them- 
selves, making SF-IETS a tool for preparing, manipulat- 
ing and reading a quantum spin in the solid state. 
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